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Abstract 

There are different approaches for the hadronic B decay calculations, recently. In this paper, 
we upgrade three of them, namely factorization, QCD factorization and the perturbative QCD 
approach based on kx factorization, by using new parameters and full wave functions. Although 
they get similar results for many of the branching ratios, the direct CP asymmetries predicted by 
them are different, which can be tested by recent experimental measurements of B factories. 
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I. INTRODUCTION 



The hadronic B decays are important for testing the standard model (SM), and also for 
uncovering the signal of new physics. Understanding non-leptonic B meson decays is crucial 
for the CKM matrix elements measurements and CP violation detection. Especially the 
direct CP violation is expected to be measured in hadronic B decays. Recently, both Belle 
and BaBar claimed to find direct CP asymmetry in B° — > 7r + 7r~ and B° — > n + K~ decays 
Q]. It is the start of direct CP measurement in B physics. 

In the theoretical side, there are much more complication in direct CP asymmetry pre- 
dictions. The direct CP violation in hadronic decays requires at least two decay amplitudes 
with different weak phase and strong phase. In the standard model, the weak phase comes 
from the so called CKM matrix, which is well defined. Although the weak phase from SM is 
clean, the strong phase requires hadronic matrix element calculation, which is usually model 
dependent. 

The simplest case is two-body hadronic B meson decays, for which Bauer, Stech and 
Wirbel (BSW) proposed the naive factorization assumption (FA) in their pioneering work 
2]. Considerable progress, including generalized FA and QCD-improved FA (QCDF) 
5|, has been made. On the other hand, technique to analyze hard exclusive hadronic 
scattering was developed by Brodsky and Lepage based on collinear factorization theorem 
in perturbative QCD (PQCD). A modified framework based on factorization theorem was 
then given in [3, |£| , and extended to exclusive B meson decays in 

The predictions of branching ratios agree well with experiments in most cases, thus it is 
difficult to tell from experiments that which method is better than others. However, the 
strong phase, which is important for the CP violation prediction, is quite sensitive to various 
approaches. The mechanism of this strong phase is quite different for various method, and 
give quite different results. The recent experimental results can make a test for the 
validity of these approaches. 

In this paper, we will first introduce the factorization approach in next section. The QCD 
factorization and improved PQCD approach based on kx factorization are then introduced 
in section IIHI and section IIV1 respectively. In section we upgrade the numerical results 
of these approaches using newest parameters. We compare the three major approaches to 
show the difference of direct CP asymmetry predicted by them. Finally the summary is 
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presented. 



II. NAIVE AND GENERALIZED FACTORIZATION APPROACH 

The calculation of non-leptonic decays involves the short-distance Lagrangian and the 
calculation of hadronic matrix elements which are model dependent. The short-distance 
QCD corrected Lagrangian is calculated to next-to-leading order 12]. In 1987, Bauer, 
Stech and Wirbel first calculate the hadronic B and D decays using the naive factorization 
approach j^. In their factorization method, the hadronic matrix element is expressed as a 
product of two factors {hih2[H e f f\B) = (hi\ Ji\B) (Z^l^lO). The first factor is proportional 
to the B — > hi form factor, while the second one is proportional to the decay constant of hi 
meson. All the perturbative channel dependent part is described by the Wilson coefficients 
of four quark operators. 

These effective four-quark operators are induced by the weak interaction in the quark 
level, which gives the short distance contribution to the non-leptonic decays of B mesons. 
The effective Hamiltonian for the charmless non-leptonic B decays is ^ 



Gf r / 10 

^ff = V q , b V*, g C, Oi + CgOg 



(1) 



where q = d,s and V q > q denotes the CKM factors. The operators 0\, O2 are tree level current 
operators. The operators O3, . . . , 0§ are QCD penguin operators. O7, . . . , O±o arise from 
electroweak penguin diagrams, which are suppressed by a/a s . Only O9 has a sizable value 
whose major contribution arises from the Z penguin. The Cj-'s are the Wilson coefficients 
of four quark operators with QCD corrections. 

Although this is a very simple method, later experiments show that many of the decay 
branching ratios explained well by the FA Q], especially for the color enhanced decays 
class I, III and class IV |3j. When charmless B meson decays are considered, efforts have 



been made to generalize the FA approach |3j, |4( . To explain the non-factorizable dominated 
class II and class V decays, phenomenological parameter is introduced. Most_of the 
branching ratios of hadronic B decays agree well with experiments by N%** = 2 



An effort is also made to predict the CP violation parameters in different hadronic B 



decays 
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Since direct CP violation requires a strong phase difference between decay 



amplitudes in addition to a weak phase difference, the precision of strong phase calculation 
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FIG. 1: The perturbative charm quark loop diagram generating strong phases in FA and QCDF. 



is essential. In naive FA, the hadronization is described by the form factor only, therefore, 
no strong phase is given. 

In generalized FA Q|, the next-to-leading order Wilson coefficients contain strong phases 
generated by the charm quark loop, where charm quark can be on shell. The diagram shown 
in Figgis also called BSS mechanism [15]. The size of this kind of strong phase is sensitive 
to the momentum of the gluon connecting to the charm quark loop. However, in FA, this 
momentum is not well defined, since all hadronic dynamics are defined only by form factors. 
In ref. the authors use k 2 = m 2 /2 ± 2GeV 2 to give the CP asymmetry parameters of 
many channels. 



III. QCD FACTORIZATION APPROACH 

In 1999, Beneke, Buchalla, Neubert, and Sachrajda (BBNS) proposed a formalism for 
two-body charmless B meson decays In this approach, they expand the hadronic matrix 
element by the heavy b quark mass and a s 



(7T7T 



Oi\B) = (7r|j 1 |B)(7r|j 2 |0) \l + ^2r n a: + 0(A QCD /m b )} . (2) 



Here O, is a local operator in the effective Hamiltonian and jxp are bilinear quark currents. 
By neglecting the power corrections of O(Aqcd), one need only calculate the order a s correc- 
tions including the vertex corrections for the four quark operators and the non-factorizable 
diagrams. These diagrams are shown in Fig|21 The first 6 diagrams have already been 
included in the generalized FA approach as next-to-leading order QCD corrections to local 




FIG. 2: Figures calculated in the QCD factorization approach. 



four quark operators. What new are the last two non-factorizable diagrams, which has a 
hard gluon line connecting the four quark operator and the spectator quark. 

They claimed that factorizable contributions, for example, the form factor F BlT in the 
B — > 7T7T decays, are dominated by long distance contributions. Hence, it is treated in 
the same way as FA, and expressed as products of Wilson coefficients and form factor F Bv . 
The non-factorizable contributions calculated perturbatively, are written as the convolutions 
of hard amplitudes with three (B, ir, ir) meson wave functions. Annihilation diagrams are 
neglected as in FA. Values of form factors at maximal recoil q 2 = m 2 and non-perturbative 
meson wave functions are all treated as input parameters. It is easy to see from eq.(J2J, that 
this equation is only applicable for those color enhanced decay modes, where the factorizable 
contribution dominates the final results. While for the color suppressed modes, where the 
non-factorizable contributions are not small, the expansion of eq.(j2J) is not right, since the 
large non-factorizable contribution is grouped into the next-to leading order term of eq.Q. 

The numerical results show that the theory and experiments agree well for those class I 
and IV decays, which are color enhanced and dominated by the factorizable contribution. 
This also agrees with the FA result, since the dominant part in eq.(0) is the same as the 
FA. The success of QCD factorization is that one can calculate the sub-leading 0(a s ) non- 
factorizable contribution (second term in eq.©) using perturbative QCD. While in the FA, 
one has to input a free parameter iV^" to accommodate the non-factorizable contribution. 

In the QCD factorization calculations, people found that there exist endpoint divergence 
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in the annihilation diagram calculations 16]. Logarithm divergence occurred at twist 2 
contribution, and linear divergence exists in twist 3 contribution. If not symmetric wave 
function, like meson, there is also soft divergence in the non-factorizable diagrams. It 
is very difficult to treat these singularity in the BBNS approach. A cut-off is introduced to 
regulate the divergence, thus makes the QCD factorization approach prediction parameter 
dependent especially in the annihilation diagrams. 

As for the strong phase in BBNS, like in FA, it mainly comes from the BSS mechanism. 
Here the momentum of the inner gluon is well defined, since wave functions are introduced. 
However, it predicts too small strong phase, because of the small gluon momentum. Hence, 
small direct CP asymmetry is predicted. There is also another source of strong phase from 
the annihilation diagrams, but strongly depends on the cut-off parameter. The strong phase 
in QCD factorization can be large due to this cut-off. 

In a word, the QCD factorization approach is at least one step forward from Naive Factor- 
ization approach. It gives systematic prediction of sub-leading non-factorizable contribution 
for the class I and class IV decays, which are dominated by the factorizable contribution. 
Problem remained is the endpoint singularity in higher order calculations, but may be solved 
with Sudakov form factors like PQCD approach. 



IV. FORMALISM OF PQCD APPROACH 

In this section, we will introduce the idea of PQCD approach. The three scale PQCD 
factorization theorem has been developed for non-leptonic heavy meson decays ^(|, based 
on the formalism by Brodsky and Lepage and Botts and Sterman []|. In the non- 
leptonic two body B decays, the B meson is heavy, sitting at rest. It decays into two 
light mesons with large momenta. There must be a hard gluon to kick the light spectator 
quark (with small momentum) in the B meson to form a fast moving light meson. So the 
dominant diagram in this theoretical picture is that one hard gluon from the spectator 
quark connecting with the other quarks in the four quark operator of the weak interaction. 
Unlike the usual factorization approach, the hard part of the PQCD approach consists of six 
quarks rather than four. We thus call it six-quark operators or six-quark effective theory. 
There are also infrared (soft and collinear) gluon exchanges between quarks. Summing 
over those leading soft contributions gives a Sudakov form factor, which suppresses the soft 
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contribution. Therefore, it makes the PQCD reliable in calculating the non-leptonic decays. 
With the Sudakov resummation, we can include the leading double logarithms for all loop 
diagrams, in association with the infrared contribution. 

There are three different scales in the B meson non-leptonic decays Myy , mi, and 1/b. The 
first scale describe the intrinsic electroweak decay of B meson, through charged current. The 
second scale mj denote the scale of energy release in the decay. Since the b quark decay 
scale rrib is much smaller than the electroweak scale mw, the QCD corrections to the four 
quark operators are non-negligible, which are usually summed by the renormalization group 
equation. This has already been done to the leading logarithm and next-to-leading order 
for years [3]. The third scale 1/b involved in the B meson exclusive decays is usually 
called the factorization scale, with b the conjugate variable of parton transverse momenta. 
The dynamics below 1/b scale is regarded as being completely non-perturbative, and can 
be parameterized into meson wave functions. The meson wave functions are not calculable 
in PQCD. But they are universal, channel independent. We can determine them from 
experiments, and they are constrained by QCD sum rules and Lattice QCD calculations. 
Above the scale 1/6, the physics is channel dependent. We can use perturbation theory to 
calculate channel by channel. 

With all the large logarithms resummed, the remaining finite contributions are absorbed 
into a perturbative b quark decay sub-amplitude H(t). Therefore the three scale factoriza- 
tion formula is given by the typical expression, 



where C(t) are the corresponding Wilson coefficients, $(a;) are the meson wave functions and 
the variable t denotes the largest mass scale of hard process H, that is, six-quark effective 
theory. The quark anomalous dimension 7^ = —q: s /tt describes the evolution from scale t 
to 1/b. Since logarithm corrections have been summed by renormalization group equations, 
the above factorization formula does not depend on the renormalization scale \i explicitly. 

As shown above, in the PQCD approach, we keep the kx dependence of the wave function. 
In fact, the approximation of neglecting the transverse momentum can only be done at the 
non-endpoint region, since kx <C k + is qualified at that region. At the endpoint, k + — > 0, 
kr is not small any longer, neglecting is a very bad approximation. By, keeping the kr 
dependence, there is no endpoint divergence as occurred in the QCD factorization approach, 




(3) 
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FIG. 3: Diagrams for B — > P1-P2 decay in perturbative QCD approach. The factorizable diagrams 
(a),(b), non-factorizable (c), (d), factorizable annihilation diagrams (e),(f) and non-factorizable 
annihilation diagrams (g),(h). 



while the numerical result does not change at other region. Furthermore, the Sudakov form 
factors suppress the endpoint region of the wave functions. 

The main input parameters in PQCD are the meson wave functions. It is not a surprise 
that the final results are sensitive to the meson wave functions. Fortunately, there are 
many channels involve the same meson, and the meson wave functions should be process 
independent. In all the calculations of PQCD approach, we follow the rule, and we find that 
they can explain most of the measured branching ratios of B decays. For example: B — > tttt 
decays, B — > np, B — ► ttuj decays [l^], B — > Kir decays [3], B — > KK decays the form 
factor calculations of B —> n, B — > p j^J, B — > Krj^ decays J21]], B — > K(p decays |3| etc. 

We emphasize that non-factorizable (FigOfc)(d)) and annihilation diagrams (FigOfe-h)) 
are indeed sub-leading in the PQCD formalism as Mb — ► 00. This can be easily observed 
from the hard functions in appendices of ref . ^[ Q] . When Mb increases, the B meson wave 
function enhances contributions to factorizable diagrams. However, annihilation amplitudes, 
being independent of B meson wave function, are relatively insensitive to the variation of 
Mb- Hence, factorizable contributions become dominant and annihilation contributions are 
sub-leading in the Mb — > 00 limit (2^. Although the non-factorizable and annihilation 
diagrams are sub-leading for the branching ratio in color enhanced decays, they provide the 
main source of strong phase, by inner quark or gluon on mass shell. The BSS mechanism 
should also be present in the PQCD approach. However this mechanism of strong phase is 
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TABLE I: Direct CP asymmetries calculated in FA [Ij, QCDF [5j and PQCD 



Q, QCDF fl and PQCD [l? 



181 for B -> -kit 



and B — > Kit decays together with the averaged experimental results at percentage. 



Quantities 


FA 


QCDF 


PQCD 


Data 


B° -» vr+vr- 


-5 ±3 


-6 + 12 


+30 ± 20 


+37 + 11 


B° -» ir+K- 


+10 ±3 


+5 + 9 


-17 + 5 


-10.9 + 1.9 


B + _> K 7T+ 


+1.7 ±0.1 


+1 + 1 


-1.0 + 0.5 


-2.0 + 3.4 


B + K+^o 


+8 + 2 


+7 + 9 


-13 + 4 


+4 + 4 



small, since it is at next-to-leading order 0(a s ) corrections. 



V. NUMERICAL RESULTS AND DISCUSSION 



In numerical analysis, we use recently updated decay constants and form factors for 
the FA and QCDF approaches. The branching ratios changes a little bit, since they are 
sensitive to these parameters. But the direct CP asymmetry changes very little. In the 
PQCD approach, we use the full set of light meson wave functions, including two twist 3 



17 



. The threshold 
[23]. The 



distribution amplitudes, where only one used in the previous papers 
resummation for the endpoint of the hard part calculation is also newly included 
numerical results for the branching ratios and CP asymmetries change only a little bit. 

As discussed in the previous section, the strong phase generated from PQCD approach 
is quite different from the FA and QCDF approaches. The direct CP asymmetry in SM is 
proportional to the sine of the strong phase difference of two amplitudes. Therefore the direct 
CP asymmetry will be different if strong phase is different. The predicted CP asymmetry 
by the three methods are shown in table HJ It is easy to see that the FA |lj] and QCDF p| 
results are quite close to each other, since the mechanism of strong phase is the same for 
them. 



Recently the two B factories measure some channels with non-zero direct CP asymmetry 



which are shown in tableffl[24]. It is claimed that direct CP has been found in B 



IX 1 7T 



and B° — > tt + K~ deca ys w ith more than 4<r signal. It is easy to see that our PQCD results of 
direct CP asymmetry [l2( agree with the experiments, especially for the experimentally well 
measured channels B° — > 7t + tt~ and B° — > n + K~ decays. Although FA and QCDF are not 
yet ruled out by experiments, but the experiments at least tell us that the dominant strong 
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phase should come from the mechanism of PQCD not QCDF. Charm quark loop mechanism, 
which gives the central value of strong phase in QCDF, is argued next-to leading order in 
PQCD. This argument is now proved by B factories experiments. 

In summary, branching ratios of both FA and QCDF depend on the values of form 
factors. The strong phase generated from these approaches is dominantly from perturbative 
BSS mechanism. In the PQCD approach, dominant by short distance contribution, the 
form factors are sensitive to meson wave functions. By including the fcr dependence and 
Sudakov suppression, there is no endpoint divergence. In the PQCD formalism, annihilation 
amplitudes are of the same order as non-factorizable ones in powers of 1/Mb, which are both 
0(1/(MbAqcd))- The strong phase comes mainly from the annihilation and non-factorizable 
diagrams in PQCD approach, which is quite different from the FA and QCDF approaches. 
The experimentally measured direct CP asymmetry implies that PQCD gives at least the 
dominant strong phase than other approaches. This will be further tested by experiments. 
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